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ABSTRACT 



We show that the sharp cutoff in the hard X-ray spectrum of NGC 4151, unusual for Seyfert 1 galaxies, can be 
reconciled with the average Seyfert 1 spectrum if we assume that the central source is completely hidden from 
our line of sight by the thick part of the accretion disk or by the broad emission-line clouds. The observed X-ray 
radiation is produced by scattering of the Seyfert 1 type spectrum in the higher, cooler parts of the accretion disk 
corona, or in a wind. A sharp cutoff appears as a result of the Compton recoil effect. This model naturally explains 
a discrepancy regarding the inclination of the central source, inferred to be low (face-on) from observations of 
the iron Ka emission line, but inferred to be high on the basis of optical and UV observations. 

Subject headings: accretion, accretion disks — galaxies: individual (NGC 4151) — galaxies: Seyfert — 
gamma rays: theory — radiation mechanisms: thermal — X-rays: galaxies 


1. INTRODUCTION AND CONCLUSIONS 

The brightest Seyfert galaxy in X-rays, NGC 4151, has a 
significantly different spectrum from the average Seyfert 
I spectrum. The average hard X-ray Seyfert 1 spectrum is well 
described by a power law with exponential cutoff at energy, 
E ( ~ 300-1000 keV, and a Compton reflection component. 
The spectrum of NGC 4151 has a much sharper decline in 
hard X-rays and no clear signature of a reflection component. 
While the X-ray spectra of both an average Seyfert 1 and NGC 
4151 can be well described by models invoking thermal 
Comptonization of the soft radiation from the accretion disk 
in a hot corona, the corona in NGC 4151 is required to be 
much thicker (Thomson optical depth r T ~~ 2) and much 
cooler ( T ( . ~ 40-50 keV) than the corona of Seyfert Is, for 
which r, - 0. 2-0.3 and T e ~ 200-300 keV (Zdziarski et al. 
1995, 1996). 

In this Letter we argue that the intrinsic spectrum of NGC 
4151 does not differ from that of Seyfert Is, if we assume that 
the direct component is hidden from our line of sight by the 
outer parts of the accretion disk or by the broad emission line 
region (BLR) close to the central source (Jourdain & Roques 
1995) and that the observed X-ray radiation is due to scattering in 
the higher, cooler parts of the accretion disk corona or in a wind. 
The observed Compton-scattered component has a much sharper 
cutoff than the intrinsic spectrum due to the Compton recoil 
effect. We show also that the primary X-ray spectrum of NGC 
4151 is consistent with thermal Comptonization in active 
regions in the vicinity of a relatively cold accretion disk and 
that optical depths and temperatures of the hot plasma do not 
differ from those of Seyfert Is. Nonthermal models cannot be 
ruled out, as a strong annihilation line will be smeared out by 
scattering. The scattered component is further filtered through 
a complex absorber. The absorption clearly visible in the X-ray 
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spectrum of NGC 4151 can be provided by the extended 
“atmosphere” of the accretion disk or BLR. 

The edge-on orientation of the NGC 4151 nucleus is 
strongly supported by the biconical geometry of the [O m] 
A5007 region (Evans et al. 1993; Pedlar et al. 1993). The 
observed geometry requires the observer to be located outside 
the cone of UV radiation that photoionizes the oxygen. Recent 
observations of the profile of the iron Ka line, showing an 
extended luminous red wing and a sharp cutoff on the blue 
side, suggest an accretion disk viewed face-on (Yaqoob et al. 
1995). This can be reconciled with the edge-on geometry 
deduced from the [O in] A5007 image, if the central source is 
observed through the radiation scattered by electrons in an 
extended corona or wind. 

Finally, after correcting the UV and X-ray luminosities for 
dilution due to scattering, one also finds that NGC 4151 has 
luminosity ratios, L i:v /Loiii and L x /L ouu typical of Seyfert Is 
(Mulchaey et al. 1994). 

2. SCATTERING MODEL 

The scattering region is assumed to be situated along the 
axis of the accretion disk in the form of a cone (we call it the 
“scattering cone”). In order to avoid additional parameters, we 
assume that the radiation coming from the central source is 
scattered to our line of sight at a given angle. This angle was 
fixed at i = 65°, which is the best estimate obtained from 
optical and radio observations (Evans et al. 1993; Pedlar el al. 
1993). 

For an arbitrary intrinsic radiation spectrum described by 
the photon number density, f nlr (x), the Compton scattered 
component can be found using the Klein-Nishina formula: 

x T S1 .[1 + (LT + XV, (1 - (1) 

where x = hv/m,.c 2 is the dimensionless photon energy, 
/x = cos /, jcj = jc/[ 1 — jt( 1 — /x)], and t* is the Thomson 
optical depth of the scattering material. Here we assume that 
the electron temperature is much smaller than m r c 2 and that 
the bulk velocity in the scattering region is small compared 
with the speed of light. Under such assumptions the scatterer 
can be treated as static and the Compton scattering cross 
section for cold electrons can be applied. Then, for small 
photon energies (.v S 0.1) the shape of the spectrum remains 
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TABLE 1 


Bfst-Fit Paramltfrs for 1991 Junf;/Jui.y (Obs 1) and 1993 May (Obs 2) 
Obshrvations of NGC 4151 
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Nortv — CPLR: exponentially cutoff power-law plus Compton reflection. Reflection model of Magdziarz 
& Zdziarski 1995 was used. PBC: pill-box-corona model by Poutanen & Svensson 19%. All errors arc for 
A* 2 =2.7. 


the same, but at higher energies, the cutoff appears much 
sharper than in the intrinsic spectrum due to the Compton 
recoil. 

3. SPECTRAL FITTING 

We fitted nearly simultaneous observations of NGC 4151 by 
ROSAT, Ginga , and OSSE in 1991 June/July (Obs 1), and 
ASCA and OSSE in 1993 May (Obs 2). Description of the 
broadband data can be found in Zdziarski et ah (1996). We 
used XSPEC version 8.5 (Shafer, Haberl, & Arnaud 1991) to 
fit the data. 

3.1. Exponentially Cut-off Power Law 
with Compton Reflection 

We apply the scattering model described in § 2, where the 
intrinsic spectrum is represented as a sum of an exponentially 
cut-off power law and a Compton-reflected component. This 
type of spectrum gives a good description of the spectra of 
Seyfert 1 galaxies (Nandra & Pounds 1994; Zdziarski et al. 
1995). We use the Compton reflection model of Magdziarz & 
Zdziarski (1995). The parameters of the model are spectral 
energy index, a, cut-off energy, E t , and the amount of 
reflection, R. For an isotropic source illuminating a flat 
reflecting slab, R = 1 is expected. The iron line component is 
modeled by the “disk-line” model of Fabian et al. (1989). The 
inner and outer radii of the accretion disk, r ( and r„, are fixed 
at 6 and 1000 Schwarzschild radii, r G , respectively. The line 
energy, £ Fc , and the emissivity index, q (characterizing the line 
emissivity as oc r“ £/ ), are allowed to vary. The intrinsic spectrum 
corresponds to the inclination / = 0°. The scattered compo- 
nent is further transmitted through a complex absorber, which 
we approximate by a conventional dual absorber model with 
column densities covering the whole source, and N 2 U 
covering a fraction C F of the source (Weaver et al. 1994). In 
our case, by the source we mean the scattered component. We 
use the abundances from Morrison & McCammon (1983). As 
discussed in Warwick, Done, & Smith (1995) and Zdziarski et 
al. (1996), the data below 1 keV show a separate soft compo- 
nent, which we model by a power law with exponential cutoff 
at 100 keV. An additional neutral absorber with column 


density, , exceeding the Galactic value of 2.1 X 10 2 ° cm 2 
(Stark et al. 1992) covers the whole source. 

In spectral fitting we use only one reflection component with 
i = 0°, which corresponds to the radiation reflected from the 
accretion disk and scattered in the scattering cone. However, 
we expect the existence of another component arising, e.g., 
from the far wall of the obscuring medium if the latter has the 
form of an opaque torus (Ghisellini, Haardt, & Matt 1994; 
Krolik, Madau, & Zycki). It is impossible to distinguish 
between these reflection components, and it is also quite 
difficult to separate reflected and scattered radiation at high 
energies, as both have a sharp cutoff due to Compton recoil. 

The Ginga data (Obs 1) provide good constraints on the 
amount of reflection, giving a best-fit value of R — 
0.30 ± 0.13, which is less than is typically found in Seyfert Is. 
In § 4 we give possible explanations for the weak reflection in 
NGC 4151. The ASCA data (Obs 2) do not constrain the 
amount of reflection, and the inferred values of R , a, and E r 
are strongly correlated. The data do not require a reflection 
component, but its presence cannot be ruled out. For Obs 2, 
we fix R at the best-fit value for Obs 1. The best-fit parameters 
are given in Table 1. 

The cutoff of the intrinsic spectrum (E c ~ 200-500 keV) 
appears at a much higher energy than in the models without 
scattering ( E { ~ 80-140 keV; Zdziarski et al. 1996), and is in 
very good agreement with the cutoff of the average Seyfert 
1 spectrum (E r = 560*2™ keV) and the cutoff of the average 
spectrum of all Seyfert galaxies (E c = 320 + jo' keV, see Zdzi- 
arski et al. 1995). Both observed spectra of NGC 4151 are 
harder than the average Seyfert 1 spectrum. For Obs 1, the 
power-law index, a, of 0.7 (cf. Table 1) is inside the range of 
observed power-law indices of Seyfert Is (Nandra & Pounds 
1994), while a ~ 0.5 in Obs 2 is much flatter. 

3.2. Two-Phase Disk-Corona Model 

As was already mentioned, spectra of Seyfert galaxies can be 
well described by thermal Comptonization models (Haardt & 
Maraschi 1993; Stern et al. 1995). In a two-phase model for an 
accretion disk-corona, soft (blackbody) radiation from the 
accretion disk gets Comptonized by hot thermal electron 
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(-positron) coronal plasma. We consider two geometries of the 
corona: plane-parallel slab and localized active region in the 
form of a pill-box (cylinder, with height-to-radius ratio equal 
to unity). The method of Poutanen & Svensson (1996) to 
compute Comptonized spectra by taking reflection from the 
cold disk into account is applied. The spectrum scattered by 
cold matter in the outer part of the corona is found from 
equation (1), where an intrinsic spectrum is the face-on 
spectrum of the accretion disk-corona system. The parameters 
of the model arc temperature of the cold disk, T hh , coronal 
electron temperature, 7;., and the optical depth of the corona, 
r r . We fix 7h b — 10 eV, because this parameter is purely 
constrained by the data. The same models for the iron line, 
absorber, and soft component as in § 3.1 are used. 

For the slab corona, the best-fit parameters (7^., 
ti) = (210 keV, 0.24) and (218 keV, 0.30) for Ohs 1 and 2, 
respectively, do not satisfy the energy balance between hot and 
cold phases (Haardt & Maraschi 1993; Stern et al. 1995). The 
pill-box corona model gives T ( = 210 keV and r, = 0.54 for 
Obs I (see Table 1 ). This is consistent with the situation where 
all energy is dissipated in an active region detached from the 
cold disk at approximately a height equal to half of its radius. 
In this condition, the X-ray source is photon starved, and the 
covering factor (fraction of the reprocessed radiation return- 
ing to the active region) is M).4 (see Stern ct al. 1995; 
Zdziarski et al. 1996). The best fit to Obs 2 gives T ( . = 290 keV 
and T; = 0.49. The energy balance requires dissipation of 
energy in an active region detached from the cold disk at a 
height comparable to its radius, implying a covering factor 
~0. 15. For pure pair corona, the compactness parameter, / 
(for definition, see, e.g., Stern et al. 1995), is about 100 in both 
cases. The corresponding model spectra arc shown in Figure 1. 

4. DISC USSION 

Observations of optical emission lines from the narrow-line 
region in NGC 4151 indicate that the gas there is illuminated 
by an ionizing continuum stronger than the continuum ob- 
served from Earth by a factor 13 (Penston et al. 1990). NGC 
4151 also has relatively small luminosity ratios L liV IL {m and 
L x /L om , compared to Seyfert 1 galaxies (Mulchaey ct al. 
1994). This suggests that intrinsic UV and X-ray luminosities 
are underestimated by at least a factor of 10. Observations of 
the “ionizing cone” (Evans et al. 1993) clearly show the 
anisotropy of UV radiation in NGC 4151. Our models, in 
which the central source is completely covered by optically 
thick matter and the only radiation we see is the radiation 
scattered in the cone, easily explain these properties. 

Let us assume that the inner edge of the scattering cone is 
r u = 30 r (l = 3 X 10 M cm (for a black hole mass 5 x 1() 7 M. ; 
Ulrich et al. 1984; Clavel et al. 1987). This is supported by the 
X-ray variability timescale (Yaqoob & Warwick 1991). The 
scatterer can be produced by a wind from the accretion disk. It 
can be highly ionized matter of normal composition or pure 
pair plasma. From observations we cannot distinguish between 
these alternatives. The ionization parameter, £ = L/r 2 n ( . , is of 
order 10 at r 10 1 " cm. Thus, highly ionized matter can 
scatter radiation out from the scattering cone without imprint- 
ing a spectral line signature on it. Assuming that the electron 
density decreases along the cone as n c = n[ ] (r/r u ) 1 and that 
r sc ~ 0.2, we can estimate the electron density at the inner 
edge of the cone, n " = 1(T cm \ The velocity of the wind 
would be approximately v ~ 0.1c, which is the terminal veloc- 


ity of a particle if radiation and gravitational forces are of the 
same order and the particle is at rest at 10 r ( , from the center. 
The rate of pair production to feed the wind, « p;iir , should be 
about 10 4S s '. Assuming L x = 3 X 1() 44 ergs s 1 and a pair 
yield of 1% (the fraction of luminosity converted into rest mass 
of pairs; see, e.g., Maciotek-Niedzwiecki, Zdziarski, & Coppi 
1995), we get a pair production rate h pitil . ~ 2 X 10 4S s \ This 
shows that a pair wind is consistent with both the inferred 
intrinsic X-ray flux and the required scattering optical depth. 
Pairs are expected to have a Compton temperature of order 
10 7 -l()‘ s K. 

One of the most intriguing properties of NGC 4151 is that 
different column densities of absorbing material are inferred 
from observations in UV (jV m ~ 10 |y -!0 21 C m Kriss et al. 

1992) and in X-rays (N u ~ 1() :2 -1()^ cm 2 ; Yaqoob ct al. 

1993) . Coexistence of warm and cold absorbers can explain 
this discrepancy (Evans et al. 1993; Warwick et al. 1995). The 
observed UV and X-rays can be scattered at similar distances 
from the central source, consistent with the observed temporal 
correlations between the flux in these two frequency bands 
(Perola et al. 1986; Edelson et al. 1996). 

The Ginga data show that the reflection component is rather 
weak in NGC 4151 (see §3.1 and Maisack & Yaqoob 1991). 
Gondek et al. (1996) found that the average Seyfert 1 also 
exhibits a deficit of reflection, having R = 0.67 ’jj;}?. If the 
underlying continuum is produced by Comptonization, then 
the spectral break appears at energies corresponding to the 
maximum of the second scattering order and is caused by 
anisotropy of the soft photons coming from the accretion disk 
(Stern et al. 1995). For a source compactness / ~ 100, the 
anisotropy break lies in the 2-10 keV energy range (see Fig. 3 
in Stern et al. 1995, and dotted curves in our Fig. 1). The 
overall spectrum, being the sum of the Compton reflected 
spectrum and a broken power law, is almost a perfect power 
law in the 2-20 keV band. This explains a deficit of reflection 
in Seyfert Is. If the intrinsic face-on spectrum has an amount 
of reflection (from the accretion disk) R ~ 0.6-0. 7, then the 
scattered radiation is expected to have R ~ 0.3- 0.4 due to the 
angular average. 

In the case of a torus-like geometry of obscuring matter, we 
can expect also some contribution due to reflection from the 
side of the torus opposite to the observer or from the outer 
part of the accretion disk. The actual fraction of the reflected 
radiation depends on the angular distribution of the intrinsic 
radiation, the exact geometry of the absorbing (reflecting) 
medium, dumpiness of the matter, etc. It can be about 15% 
for a torus geometry, assuming that optically thick material 
just covers the central source and the inclination angle is 65°, 
and less than 4% if the surface of the accretion disk is conelike 
with a constant opening angle ~120°. This would correspond 
to the amount of reflection R ~ 0. 1-0.4 assuming r*. ~ 0.3 
(note that in our model the intrinsic spectrum corresponds to 
/ = 0°, and R is normalized to the scattered component). 
Since the data suggest R < 0.5, we conclude that either 
geometrical effects reduce the amount of reflection from the 
torus, the central source is anisotropic, and/or the reflecting 
(obscuring) matter is clumpy. 

The equivalent width of the iron line (EW ~ 200-300 eV) is 
in agreement with the predictions of the two-phase disk- 
corona models (Poutanen, Nagendra, & Svensson 1996). It is 
significantly higher than expected from the irradiation of cold 
matter by isotropic X-ray radiation having a power-law shape 
(George & Fabian 1991), due to the anisotropy of the Comp- 
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Fig. 1. — (a) Spectrum of NGC 4151 observed in 1991 June/July by ROSAT , 
Ginga, and OSSE. Solid curve represents the best-fit model spectrum for 
pill-box-corona model (PBC; see Table 1 ). Dashed curve represents the intrinsic 
spectrum of the disk-corona system which is seen by the scattering medium 
(normalized to the scattered component). Dotted curve represents the under- 
lying Comptonized spectrum (without reflection from the cold disk), which can 
be approximated by a broken power law. The upper limits are 2 a. {h) Same as 
(a ), but for observations in 1993 May by ASCA and OSSE. 
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tonized radiation. While the inclination of the central torus- 
like region is suggested to be high based on the Balmer line 
reverberation mapping (Maoz et al. 1991), recent observations 
of the iron line profile (Yaqoob et al. 1995) strongly suggest 
that we see the accretion disk almost face-on. This discrepancy 
has a simple explanation in the scattering model, as we see the 
central source via a “mirror” situated above the disk. 

The intrinsic spectrum of NGC 4151 can be produced by 
thermal Comptonization in localized active regions above a 
cold accretion disk. The compactness parameter is high, and 
therefore a significant fraction of hot plasma is expected to be 
in the form of electron-positron pairs. A pair model can also 
naturally explain the origin of the scattering material, as a pair 
wind. Despite the fact that a thermal model gives an excellent 
fit to observations, the presence of nonthermal processes 
cannot be ruled out. While nonthermal models predict a large 
flux at E ~ 500 keV due to pair annihilation, this flux is 
significantly diminished due to scattering by cold matter. 

An increasing fraction of the central source radiation is 
transmitted through the obscuring matter due to the Klein- 
Nishina decline of the Compton scattering cross section at 
higher energies. Upper limits on the y-ray flux at E ~ 400 keV 
(see Fig. 1) provide some constraints on the Thomson optical 
thickness of the obscuring matter, r obsc . Assuming that the 
intrinsic flux is 10 times larger than the observed one, we get 

Tthsc ~ b. 

A crucial test for the scattering model will be future 
observations of polarization in X-rays with the Spectrum-X-y 
satellite. The scattering model predicts a polarization of about 
20%-30% perpendicular to the axis of the scattering cone. 
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